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ABSTRACT  1 

We  examine  the  imaging  of  standard  Brownian  Fractal  surfaces,  and 
find  that,  given  certain  assumptions,  a  Fractal  surface  with  power  spectrum 
proportional  to  f~@  has  an  image  with  power  spectrum  proportional  to 

f2~P. 


1  Introduction 

Fractal  Brownian  functions  have  found  widespread  use  in  computer  graphics 
because  they  produce  surfaces  that  closely  resemble  natural  surfaces.  It  is 
natural,  therefore,  to  ask  if  we  may  apply  such  an  attractive  model  of  nat¬ 
ural  surface  shape  to  the  problems  of  image  analysis.  Pentland  [l]  showed, 
by  modeling  natural  surfaces  using  Brownian  Fractals,  that  interesting  in¬ 
formation  about  the  imaged  surface  could  be  obtained. 

In  his  analysis  Pentland  modeled  the  imaged  surface  using  a  vector¬ 
valued  plane-to-surface-normal  Brownian  Fractal  function.  The  non-constructive 
nature  of  this  model,  however,  places  unfortunate  restrictions  upon  its  use¬ 
fulness,  In  this  paper  we  analyze  the  imaging  of  surfaces  modelled  by  Frac¬ 
tal  Brownian  plane-to-elevation  functions  of  the  sort  standardly  used  in 
computer  graphics.  We  thus  obtain  a  model  of  surface  and  image  that  is 

'This  research  was  made  possible  by  National  Science  Foundation,  Grant  No.  DCR- 
85-19283,  by  Defense  Advanced  Research  Projects  Agency  contract  no.  MDA  903-83-C- 
0027,  and  by  grants  from  the  Systems  Development  Foundation  and  the  Alfred  P.  Sloan 
Foundation.  We  would  like  to  thank  Gene  Switkes  and  Eugene  Wong  of  Berkeley  and 
Andy  Hanson  of  SRI  for  helpful  discussions. 
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constructive  in  nature,  and  therefore  of  potentially  greater  use  in  image 
analysis. 


2  Standard  IVactal  Brownian  Models 

In  one  topological  dimension,  a  Fractal  Brownian  line  function  V)j(t)  has 
zero-mean  Gaussian  increments  with  variance 

<[VIi(t+S)-VH{t)}*>  oc  |5|2"  (1) 

where  He( 0,  l).  H  is  related  to  the  fractal  dimension  D  of  Vji(t)  by 

D  —  2  —  H  (2) 

Further,  V//(t)  has  a  random-phase  Fourier  spectrum  with  power  F h{I) 
such  that 

F*(/)  «  rp  (3) 

H  is  related  to  /?  by 

P=2H+1  (4) 

In  two  topological  dimensions,  any  “slice”  through  an  isotropic  Fractal 
Brownian  surface  Vjj(x ,y)  is  a  Fractal  Brownian  line  function  of  identical 
H .  Thus  we  want 

<  [Vf/(3+  5  cos-/,  y+  5sin7)  —  Vj/(x,  y)]2  >  oc  |5|2i/  (5) 

independent  of  7.  Voss  [2]  shows  that  this  requires  the  two-dimensional 
power  spectrum  Fj/(/,  9)  of  the  surface  to  be 

Fid/,*)  r'3  (e) 

where 

P  =  2H  +  2  (7) 

These  two-dimensional  surfaces  can  be  constructed  by  filtering  Gaussian 
noise  [2],  or  (more  efficiently  and  flexibly)  by  constructing  the  an  appropriate 
Laplacian  Pyramid  [3].  Those  in  computer  graphics  often  use  a  recursive 
subdivision  algorithm  to  approximate  true  Fractal  Brownian  surfaces  [4]. 
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3  Imaging  of  Fractal  Brownian  Surfaces 


Let  z  —  y)  be  a  Fractal  Brownian  surface,  and  let  us  assume  that: 

(1)  the  surface  is  Lambertian  (we  will  later  partially  relax  this  assump¬ 
tion)  , 

(2)  the  surface  is  illuminated  by  (possibly  several)  distant  point  sources, 

(3)  the  surface  is  not  self-shadowing. 

We  will  also  take  z  <  0  within  the  region  of  interest,  and  assume  ortho¬ 
graphic  projection  onto  the  x,  y  plane. 

We  will  let  L  =  (cost sin <7, sinrsincr, cosct)  be  the  unit  vector  in  the 
mean  illuminant  direction,  where  r  is  the  tilt  of  the  illuminant  (the  angle 
the  image  plane  component  of  the  illuminant  vector  makes  with  the  :r-axis) 
and  a  is  its  slant  (the  angle  the  illuminant  vector  makes  with  the  z-axis). 

Then  the  normalized  image  intensity  I{x,y)  will  be 


where 


p  cos  t  sin  a  4-  q  sin  t  sin  a  +  cos  a 

(8) 

(p2  +  <72  +  l)1/2 

p=fxV„(x,y) 

(9) 

(10) 

[Note:  as  true  Fractal  Brownian  surfaces  are,  strictly,  nowhere  differen¬ 
tiable,  we  will  assume  throughout  an  approximation  smoothed  sufficiently 
to  allow  p  and  q  to  exist.  See  [5]  for  justification  of  this  approach.] 


3.1  Spectral  Properties 

If  we  then  take  the  Taylor  series  expansion  of  I(x,  y)  about  p,  q  =  0  through 
the  quadratic  terms,  we  obtain 

r/\  ...  COS  O  .  9  /  \ 

l{x,  y)  «  coser  +  pcosr  smer  +  q  sinr  sine - - — (p  +?)  (11) 

This  expression  gives  an  excellent  approximation  if  p,q  «1.  We  note 
that  for  real  surfaces,  such  as  mountains,  the  maximum  surface  slope  rarely 
is  more  than  15°,  i.e.,  typically  p2  +  q2  <  0.1.  Under  these  conditions  the 
linear  terms  of  Eqn.  (ll)  will  dominate  the  power  spectrum  except  when 
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the  average  illuminant  is  within  ±6°  of  the  viewer’s  position,  i.e.,  when 

sin <r  <  0.1. 

The  complex  Fourier  spectrum  Fjj(f,6)  of  Vjj(x>  y)  is,  from  Eqn.  (6), 

FH(f,0)  =  r0l2ei^  (12) 

where  <f>  is  a  random  variable  uniformly  distributed  on  (0,27t),  and  </>/,«  is 
the  value  “drawn”  at  position  (/,  6 )  in  the  Fourier  plane. 

Now  since  p  and  q  are  partial  derivatives  of  Vjj,  their  transforms  Fp  and 
Fq  are  related  to  Fjj  in  an  elementary  fashion.  We  can  write 


Fp(f,6)  =  2^  cos  0/1-^/2e«'W/.«+^2) 

(13) 

Fg{f,0)  =  2v  sin  (*/,<! +V2) 

(14) 

3.1.1  Case  1 

When  p,q  are  small  and  the  illuminant  is  not  behind  the  viewer  (e.g.,  sine  > 
0.1)  then  we  may  neglect  the  quadratic  terms  of  Eqn.  (11)  and  consider 

Ii  (x,y)  =  cos  or  +  p  cos  t  sin  or  +  q  sin  r  sin  or  (15) 

In  this  case,  the  Fourier  transform  of  the  image  I\  is  (ignoring  the  DC  term): 

Fi1(f,6)  =  2tt  sinrr/^V^+^fcos  0  cos  r  +  sin  9  sin  t]  (16) 

and  the  power  spectrum  is 

F (/,  9)  =  4tt2  sin2  ct/2~^[cos0  cosr  +  sin  9  sin  r]2  (17) 

This  spectrum  depends,  as  expected,  upon  the  illuminant  direction.  As 
with  the  Fractal  surface  itself,  however,  the  spectral  falloff  is  isotropic:  the 
log  of  the  power  spectrum  of  the  image  has  slope  2-/3  with  respect  to  log 
frequency  at  almost  all  orientations  (excepting  a  set  of  measure  zero  where 
9  =  t  ±  jt/2). 

3.2  Case  2 

When  the  mean  illuminant  vector  is  almost  parallel  to  the  viewing  direction 
(i.g.,  sin  o  w  0)  the  quadratic  terms  of  Eqn.  (11)  can  dominate  and  the 
image  of  a  fractional  Brownian  surface  will  look  like 

h{x,y)  =  cos<r(l-  (p2  +  q2)/2)  (18) 
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To  within  a  constant  factor,  and  ignoring  DC,  the  power  spectrum  Fj2  of  this 
image  will  be  the  Fourier  transform  of  the  autocorrelation  function  J?p2+?2 
of  p2  +  <72,  which  we  now  investigate.  The  autocorrelation  function  is  written 


^p2+(J=(z,y)  =  <  [?(**,  I/*)2  +  g(z*,y*)2] 

[p(x*  +  x,  y*  +  y)2  +  q{x  +  x,  y*  +  y)2]  > 

-  <p(x*,y*)2  +  g(x*,y*)2  > 

<  p(x *  +  x,  y*  +  y)2  +  ?(**  +  x,  y*  +  y)2  ^19) 


where  <>  denotes  expectation.  Since  p  and  q  are  stationary,  we  can  write 

^p2+92  (-c>  y) 

=  <  PooPxy  >  +  <  qlo<lly  > 

+  <  9ooPxy  >  +  <  Poorly  > 

-  <  p200  ><  p\y  >  -  <  9q0  ><  q2xy  > 

-  <  ?oo  ><  ply  >  -  <Poo><  ?*v  > 

(20) 

where  we  have  put  pxy  for  p(x,  y),  etc.  Since  p  and  q  are  Gaussian  processes, 
the  first  four  terms  on  the  right  hand  side  of  Eqn.  (20)  can  be  rewritten 

A 

O  to 

o 

H  ^ 

V 

=  2<p0oPxv>2  +  <Poo><Pzy> 

(21) 

A 

> 

O 

o 

V 

=  2  <  <7oo?zy  >2  +  <  qoo  ><  qly  > 

(22) 

<  7ooPxy  > 

=  2  <  yooPzy  +  <  ?00  '><~  Pxy  ^ 

(23) 

A 

M  H 

Cr 

O 

V 

=  2  <  paoqXy  >2  +  <  Poo  ><  qly  > 

(24) 

Thus 

■^p!+g2  y ) 

—  2(<  poOPzy  >2  +  <  Pooqzy  > 2 

+  <  qOQPzy  >2  +  <  qooqzy  > 2) 

(25) 

We  know  what  the  Fourier  transforms  of  the  autocorrelations  of  p  and  of  q 
are;  they  are  just  the  power  spectra  of  p  and  of  q  respectively,  easily  obtained 
from  Eqns.  (13)  and  (14): 


Fp(/,  6)  =  47r2  cos2  6f2~P 

(26) 

Vq(f,6)  =  4n2sm2ef2-f3 

(27) 

To  see  what  the  Fourier  transform  of  <  pooqxy  >=<  qoo pxy  >  must 
be,  we  construct  the  expression  for  the  autocorrelation  i?j,  of  I\{x,y)  — 
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(p cos r  sin cr  4- g sin  r  sine)  and  expand  it  to  obtain: 

Rli  (*,!/)  =  cos2  t  sin2  (T  <  pooPxy  > 

+  sin2  t  sin2  er  <  qooqzy  > 

+2  cos  t  sin  t  sin2  a  <  pooqXy  >  (28) 

If  we  now  take  the  Fourier  transform  of  this  expression  we  find,  as  we  already 
know  the  power  spectrum  of  Ii  (cf.  Eqn.  (17)),  and  the  power  spectra  of  p 
and  g,  that  we  can  solve  for  the  Fourier  transform  Fpq  of  <  pooqzv  >.  We 
obtain 


Fpq  =  47r2sin0cos0/2  ^  (29) 

With  this  result,  we  find  that  the  Fourier  transform  of  I7p2+92  (x,  y)  is 
then 

Fr  =  2(Fp  *  Fp)  +  4(Fp,  *  Fpq)  +  2(F9  *  F,)  (30) 

where  *  denotes  convolution.  In  Cartesian  coordinates,  we  can  write 


OCX 

-C 


[u*(u*  -  u)  +  t;*(v*  -  u)]2 
[(u*2  +  n+2)((u*  -  u)2  +  (v*  -  u)2)]^/2 


du* dv*  (31) 


This  integral  diverges  at  u  =  v  =  0  for  any  /?.  When  /?  <  3,  it  diverges 
at  all  u, u  because  of  behavior  of  the  integrand  as  oo.  We  note, 

though,  that  the  spectral  power  of  a  real  surface  does  not  go  to  infinity  as 
frequency  goes  to  zero,  and  any  sufficiently  high  frequency  components  of  a 
surface  can  be  ignored  in  a  model  of  a  physically  realizable  imaging  process; 
thus  Fr  can  be  taken  as  well  defined  everywhere. 

To  find  an  approximation  to  Fr,  note  that  since  (3  >  2,  f2~&  becomes 
quickly  very  large  as  /  approaches  0.  Therefore  we  might  think  of  the  self¬ 
convolutions  represented  in  Eqn.  (30)  as  instead  being  convolutions  with 
a  localized  averaging  kernel,  say  a  Gaussian  with  small  variance.  Then  (at 
least  for  frequencies  sufficiently  far  from  zero) 


Fp  *  Fp  «  Fp  Fpq  *  Fpq  «  Fpq  Fq  *  Fq  «  F,  (32) 

and  thus 

logF/2  oc  log  Fr  oc  2-/3  (33) 

The  conclusion,  therefore,  is  that  when  the  mean  illumination  is  at  the 
viewer’s  position  the  image  will  have  a  power  spectrum  falloff  approximately 
proportional  to  f2~& ,  i.e.,  the  same  relationship  between  surface  Fractal 
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Figure  1:  Comparison  of  predicted  (top  row)  versus  the  actual  (bottom 
row)  radially-averaged  log  Fourier  power  versus  log  frequency  for  images 
of  Fractal  dimension  D  =  2.25,2.5,2.75,  assuming  the  illuminant  is  at  the 
viewer’s  position. 

scaling  parameter  and  image  Fractal  scaling  parameter  that  we  found  with 
oblique  illumination.  This  conclusion  is  supported  by  the  results  of  the 
following  simulation. 

3.2.1  Simulation 

Range  maps  z  =  Vu(x,y )  for  H  =  0.25,  0.5  and  0.75  (corresponding  to 
surface  of  Fractal  dimension  of  D  —  2.75,  2.5,  and  2.25)  were  illuminated 
and  viewed  along  the  z  axis  and  the  power  spectrum  exponent  was  observed 
for  the  resulting  image.  For  each  H,  Vh  was  scaled  to  have  -E(|p|)  =  0.1. 
Figure  1  shows  graphs  of  the  log  Fourier  power  versus  log  frequency  of  the 
images’  spectrum;  power  at  a  frequency  was  averaged  over  all  orientations 
at  that  frequency. 

3.3  Oblique  Viewpoints 

So  far  we  have  considered  only  views  straight  down  on  the  surface.  We 
would  like  to  consider  the  case  of  oblique  views  of  the  surface.  Without  loss 
of  generality  we  may  consider  only  the  case  where  in  the  viewer’s  co-ordinate 


7 


system  ( x*,y*,z *)  the  surface  has  partials  ( p*,q *)  such  that 


p*=p+c  q*  =  q  (34) 

where  ( p,q )  are  (as  in  Eqns.  (9)  and  (10)  the  partials  measured  in  the  sur¬ 
face’s  co-ordinate  system.  Thus  the  surface  has  an  average  slope  dz*  jdx*  =  c 
with  regard  to  the  viewer. 

In  this  case  Eqn.  (11)  becomes 


If  continue  to  assume  that  ptq  are  small,  then  in  the  oblique  co-ordinate 
system  the  frequency  /*  =  8f,  where  8  =  (1  +  c)-1  cos2  6  +  sin2  6,  and  f  is 
the  freqency  in  the  surface  co-ordinate  system  (i.e.,  there  is  foreshortening). 
Similarly,  the  orientation  6*  in  the  oblique  co-ordinate  system  is  a  function 
of  the  orientation  6  in  the  surface  co-ordinate  system.  Thus  Eqns.  (13)  and 
(14)  become: 

Fp*(f\e*)  =  27ccos6*(e)81/2-(,/if1-(,/ieittf'0+*/2')  (36) 

Fg‘(f*,6*)  =2*sin0*(0)S1/2-^V1_^/^¥+^2)  (37) 

If,  as  in  Case  1  above,  we  are  in  a  condition  in  which  we  may  ignore  the 
quadratic  terms  of  Eqn.  (35),  we  then  find  that  the  complex  spectrum  of 
the  image  (again  ignoring  the  DC  terms)  is  given  by: 


=  [cos cr  +  (p  +  c)  cost  sin tr  +  gsin  r  sin  c](l  +  c2) 


cosc  +  c cost smtr.  2 


(P2  +  ?2) 


(35) 


Fj(f*,9*)  =  2ir(l  +  cz)-1^f1-^2ei^f’B+^ 

[(1  -  c)S W-V*  cos  6* (6)  cost  +  sin  9* (6)  sin  r] 

-  cos  O*(0)  cos (77r8lf2-p/if1-p/2ei{<t,^+v/2)  (38) 

and  thus  the  power  spectrum  of  the  image,  although  now  complexly 
dependant  upon  the  illuminant  direction  and  average  angle  of  view,  still 
has  an  isotropic  falloff  proportional  to  f2~ &  as  in  the  non-oblique  case.  A 
similar  result  holds  for  the  case  in  which  the  quadratic  terms  of  Eqn.  (35) 
dominate. 
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3.4  Non-isotropic  Surfaces 

We  have  considered  only  Fractal  surfaces  with  uniform  spatial  statistics; 
however  many  natural  surfaces  (e.g.,  tree  bark)  can  be  modeled  only  with 
anisotropic  statistics.  We  would  therefore  like  to  consider  surfaces  Vfj{x,  y ) 
with  “stretched”  spectra: 

F(f,0)  =  Sf-P^j+'-o  (39) 

where  again  S  =  (1  +  c)-1cosff  +sin2  8.  Examining  Eqns.  (6)  -  (33)  using 
this  new  surface  function,  we  find  that  this  new  multiplicative  stretching  fac¬ 
tor  simply  propagates  through  our  calculations,  resulting  in  an  “stretched” 
image  spectrum  whose  spectral  falloff,  however,  is  still  isotropic  and  whose 
rate  is  the  same  as  in  the  isotropic  case. 

3.5  Occlusion,  Shadows  and  Highlights 

We  can,  to  a  limited  extent,  analyze  cases  involving  non-Lambertian  sur¬ 
faces,  self-shadowing  and  cast  shadows.  The  following  sections  are  our  first 
step  towards  a  more  complete  analysis. 

3.5.1  Occlusions  and  Shadows 

When  the  viewing  angle  becomes  sufficiently  oblique  the  high  points  of  the 
surface  will  begin  to  occlude  the  remainder  of  the  surface;  for  viewing  an¬ 
gles  that  are  extremely  oblique,  the  effects  of  such  occlusion  will  begin  to 
dominate  the  image.  Similarly,  when  the  illuminant  becomes  sufficiently 
oblique,  there  will  be  self-shadowed,  and  when  the  illuminant  direction  be¬ 
comes  nearly  tangential  to  the  surface  then  the  image  will  be  dominated  by 
the  effects  of  cast  shadows. 

When  a  surface  occludes  itself,  there  is  a  discontinuous  jump  between 
points  that  are  adjacent  in  the  image.  Because  of  this  jump,  the  orienta¬ 
tion  at  these  image-adjacent  points  will  be  only  very  weakly  correlated  thus 
producing  a  step  discontinuity  in  the  image.  The  spectrum  of  such  a  discon¬ 
tinuity  is  proportinal  to  1//;  i.e.,  the  spectrum  of  a  step  discontinuity  falls 
off  much  more  slowly  than  does  the  normal  image  spectrum.  As  a  result 
when  such  steps  are  added  to  the  image  the  high  freqency  content  of  the 
image  quickly  becomes  dominated  by  their  spectral  contribution. 

Similarly,  casts  shadows  introduce  large  amplitude  step  edges  into  the 
image.  Again,  the  high  frequency  content  of  such  edges  quickly  comes  to 
dominate  the  spectral  content  of  the  image. 
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3.5.2  Highlights 


Highlights  occur  when  the  surface  orientation  exactly  reflects  the  incident 
illumination  toward  the  viewer;  i.e.,  they  are  associated  with  a  particular 
level-sets  of  the  (independent)  p  and  q  functions.  That  is,  hightlights  occur 
at  the  intersection  of  separate  level-sets  of  the  p  and  q  functions.  From  Man¬ 
delbrot  [6],  we  know  that  both  the  level-set  of  a  Fractal  Brownian  surface 
and  the  surface  itself  have  the  same  scaling  parameter  H.  Thus  the  set  of 
image  highlights  will  be  the  intersection  of  two  sets  of  dimension  D  =  2H+1. 
We  speculate  that  the  dimension  of  this  intersection  set  is  D  =  2H.  If  so, 
then  by  identifing  the  set  of  surface  hightlights  we  can  (by  estimating  their 
Fractal  scaling  parameter)  directly  obtain  the  Fractal  scaling  parameter  of 
the  surface  itself. 

4  Summary 

We  have  analyzed  the  imaging  of  Fractal  Brownian  surfaces  such  as  are  con¬ 
structed  in  computer  graphics.  We  have  found  that,  given  the  Lambertian 
assumption  and  no  self-shadowing,  that  the  image  of  a  Fractal  Brownian 
surface  with  spectrum  proportional  to  f~&  has  a  spectrum  proportional 
to  f2~P.  Thus  the  main  result  of  this  paper  is  that,  in  the  case  of  the 
non-isotropic  Brownian  Fractals  used  in  computer  graphics,  we  can  use  the 
Fractal  scaling  parameter  (spectral  falloff)  of  the  image  to  predict  the  Frac¬ 
tal  scaling  parameter  of  the  surface.  We  also  have  shown  that  this  result 
is  not  affected  by  multiplicative  “stretching”  of  the  surface  and,  given  that 
p2  +  q 2  is  small,  is  relatively  unaffected  by  viewing  position. 

These  results  agree  with  those  obtained  using  an  isotropic  Fractal  surface 
model  [l].  The  major  advance  of  this  paper  is  that,  because  we  have  used 
a  Fractal  surface  model  that  is  defined  constructively,  we  can  now  easily 
synthesize  surfaces  whose  images  have  the  desired  Fractal  scaling  behavior. 

And  finally,  it  appears  that  self-shadowing  and  self-occlusion  produce 
image  features  which  “swamp”  the  normal  Fractal  scaling  characteristics 
with  /-1  noise.  The  fractal  scaling  parameter  of  imaged  highlights,  however, 
appears  provide  us  with  a  method  of  obtaining  an  independent  estimate  of 
the  surface’s  Fractal  scaling  parameter. 
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